Small-world networks: Evidence for a crossover picture 
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Watts and Strogatz [Nature 393, 440 (1998)] have recently introduced a model for disordered 
networks and reported that, even for very small values of the disorder p in the links, the network 
behaves as a small-world. Here, we test the hypothesis that the appearance of small-world behavior 
is not a phase-transition but a crossover phenomenon which depends both on the network size n 
and on the degree of disorder p. We propose that the average distance £ between any two vertices of 
the network is a scaling function of n/n*. The crossover size n* above which the network behaves 
as a small-world is shown to scale as n*(p < 1) ~ p~ T with r m 2/3. 



PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18.Sn 



Two limiting-case topologies have been extensively 
considered in the literature. The first is the regular lat- 
tice, or regular network, which has been the chosen topol- 
ogy of innumerable physical models such as the Ising 
model or percolation [§]-B|. The second is the random 
graph, or random network, which has been studied in 
mathematics and used in both natural and social sciences 

Erdos and co-workers studied extensively the proper- 
ties of random networks — see jijj for a review. Most of 
this work concentrated on the case in which the number 
of vertices is kept constant but the total number of links 
between vertices increases |L7|]: The Erdos- Renyi result 
p8| states that for many important quantities there is a 
percolation-like transition at a specific value of the av- 
erage number of links per vertex. In physics, random 
networks are used, for example, in studies of dynam- 
ical problems 19 2(j], spin models and thermodynam- 
ics random walks (22), and quantum chaos [p3[ . 
Random networks are also widely used in economics and 
other social sciences ||,|4],|5| to model, for example, in- 
teracting agents. 

In contrast to these two limiting topologies, empirical 
evidence |26|,^7j suggests that many biological, techno- 
logical or social networks appear to be somewhere in be- 
tween these extremes. Specifically, many real networks 
seem to share with regular networks the concept of neigh- 
borhood, which means that if vertices i and j are neigh- 
bors then they will have many common neighbors - 
which is obviously not true for a random network. On 
the other hand, studies on epidemics |l4|,[l5],^6| show that 
it can take only a few "steps" on the network to reach a 
given vertex from any other vertex. This is the foremost 
property of random networks, which is not fulfilled by 
regular networks. 

To bridge the two limiting cases, and to provide a 
model for real- world systems p8lp9[, Watts and Stro- 



gatz [ p6|]27| ] have recently introduced a new type of net- 
work which is obtained by randomizinga fraction p of the 
links of the regular network. As Ref. [|26| , we consider as 
an initial structure (p — 0) the one-dimensional regular 
network where each vertex is connected to its z nearest 
neighbors. For < p < 1, we denote these networks dis- 
ordered, and keep the name random network for the case 
p = 1. Reference |2(| reports that for a small value of 
the parameter p — which interpolates between the regu- 
lar (p — 0) and random (p = I) networks — there is an 
onset of small- world behavior. The small- world behav- 
ior is characterized by the fact that the distance between 
any two vertices is of the order of that for a random 
network and, at the same time, the concept of neigh- 
borhood is preserved, as for regular lattices [Fig.|l|. The 
effect of a change in p is extremely nonlinear as is visu- 
ally demonstrated by the difference between Figs. [j]a,d 
and Figs. |l|b,e where a very small change in the adja- 
cency matrix leads to a dramatic change in the distance 
between different pairs of vertices. 

Here, we study the origins of the small-world behav- 
ior [ p8|p9[| . In particular, we investigate if the onset of 
small- world networks is a phase transition or a crossover 
phenomena. To answer this question we consider not 
only changes in the value of p but also in the system size 
n . 

The motivation for this study is the following. In a 
regular one-dimensional network with n vertices and z 
links per vertex, the average distance t between two ver- 
tices increases as n/(2z) — the distance is defined as the 
minimum number of steps between the two vertices. The 
regular network is similar to Manhattan: Walking along 
5 th Avenue from Washington Square Park in i th Street to 
Central Park in 59 th Street, we have to go past 55 blocks. 
On the other hand, for a random network, each "block" 
brings us to a point with z new neighbors. Hence, the 
number of vertices increases with the number of steps 
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k as n ~ z , which implies that £ increase as Inn/ In z. 
The random network is then like a strange subway system 
that would directly connect different parts of Manhattan 
and enable us to go from Washington Square Park to 
Central Park in just one stop. In view of these facts, it is 
natural to enquire if the change from large- world [£ ~ n) 
to small- world (£ ~ Inn) in disordered networks occurs 
through a phase transition for some given value of p |3C| l 
or if, for any value of p, there is a crossover size n*(p) be- 
low which our network is a large- world and above which 
it is a small-world. 

In the present Letter we report that the appearance of 
the small-world behavior is not a phase-transition but a 
crossover phenomena. We propose the scaling ansatz 



£(n,p)~n*F^) , 



(1) 



where F(u <C 1) ~ u, and F(u ^> 1) <~ In it, and n* 
is a function of p |31|. Naively, we would expect that 
when the average number of rewired links, pnz/2, is 
much less than one, the network should be in the large- 
world regime. On the other hand, when pnz/2 3> 1, 
the network should be a small- world p2]. Hence, the 
crossover size should occur for n*p = 0{1), which im- 
plies n* ~ p~ T with t = 1. This result relies on the fact 
that the crossover from large to small worlds is obtained 
with only a small but finite fraction of rewired links. We 
find that the scaling ansatz (Q) is indeed verified by the 
average distance i between any two vertices of the net- 
work. We also identify the crossover size n* above which 
the network behaves as a small-world, and find that it 
scales as n* ~ p~ T with t sw 2/3, distinct from the trivial 
expectation r = 1. 

Next, we define the model and present our results. We 
start from a regular one-dimensional network with n ver- 
tices, each connected to z neighbors. We then apply the 
"rewiring" algorithm of p3 to this network. The algo- 
rithm prescribes that every link has a probability p of 
being broken and replaced by a new random link. We re- 
place the broken link by a new one connecting one of the 
original vertices to a new randomly selected vertex. Each 
of the other n — 2 vertices — we exclude the other vertex 
of the broken link — has an a priori equal probability of 
being selected but we then make sure that there are no 
duplicate links. Hence, the algorithm preserves the total 
number of links which is equal to nz/2. A quantity that 
is affected by the rewiring algorithm is the probability 
distribution of local connectivities. For p ~ 0, this prob- 
ability is narrowly peaked around z, but it gets broader 
with increasing p. For p = 1, the average and the stan- 
dard deviation of the local connectivity are of the same 
order of magnitude and equal to z. 

Once the disordered network is created, we calculate 
the distance between any two vertices of the network and 
its average value I. To calculate all the distances between 
vertices, we use the Moore-Dijkstra algorithm p3| whose 
execution time scales with network size as n 3 Tnn. We 



perform between 100 and 300 averages over realizations 
of the disorder for each pair of values of n and p. 

Here, we present results for three values of connectivity 
z = 10, 20 and 30 and system sizes up to 1000. The scal- 
ing ansatz ([!]) enables us to determine n*(j>) from £{n) 
at fixed p. Indeed, £(n 3> n*) ~ n*lnn which implies 
that n* is the asymptotic value of di/d(hm) [Fig. |(a)]. 
Figure ||(b) shows the dependence of n* on p for different 
values z. We hypothesize that 



n* - p T g(p) 

In z 



(2) 



where the term in z arises from the fact that £ = Inn/ In z 
for a random network (p — 1), and g(p — > 1) — > 0. More- 
over, g(p) approaches a constant as p — > 0, leading to 



P 



(3) 



for small p. Due to the effect of g and the fact that n < 
1000 in our numerical simulations, we are constrained to 
estimate r from the region 2.5 x 10~ 4 < p < 2 x 10~ 2 . 
For all values of z we obtain r = 0.67 ± 0.10 [Fig. §. 

Using this value of r and the scaling form (|l|) we are 
able to collapse all the values of £(n,p) onto a single curve 
[Fig. ||. This data collapse confirms our scaling ansatz 
and estimate of r. 

In summary, we have shown that the onset of small- 
world behavior is a crossover phenomena and not a 
phase transition from a large-world to a small one. The 
crossover size scales as p~ T with r ~ 2/3. The surpris- 
ing fact that t < 1, shows that the rewiring process is 
highly nonlinear and can have dramatic consequences on 
the global behavior of the network. This implies that in 
order to decrease the radius of a network it is necessary 
to rewire only a few links. We also note that the value of 
the exponent r will likely depend on the dimensionality 
of the initial regular network. This point will be adressed 
in future work. 

We believe that the disordered networks introduced in 
p6| may constitute a promising topology for more real- 
istic studies of many important problems such as flow 
in electric-power or information networks, spread of epi- 
demics, or financial systems. The results reported here 
support this hypothesis because they suggest that, for 
any given degree of disorder of the network, if the sys- 
tem is larger than the crossover size, the network will be 
in the small-world regime. 
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FIG. 1. Effect of disorder 

on the distance between vertices of the network (go to http://polymer.bu.edu/~amaral/Networks.html for color pictures). 
We consider here matrices with z = 10, n = 128 and with periodic boundary conditions, that is, vertex 1 follows vertex n. 
Adjacency matrices for (a) a regular one- dimensional network where each vertex is connected to its z nearest neighbors, (b) a 
disordered network with p — 0.01, and (c) a random network. Black indicates that a link is present between the two vertices 
while gray (green in the color picture) indicates the absence of a link. Note that (a) and (b) are nearly identical. Distance 
matrices for (d) the regular network, (e) the disordered network with p = 0.01, and (f) the random network. We use the relief 
of the surface and a color scheme to represent the distance between two vertices. Greater height indicates larger distance. The 
color scheme is the same for the relief and for the contour lines: Distance increases from very dark gray (blue) to gray (green) 
to light gray (yellow) to dark gray (red). For the regular network, the contour lines are parallel to the diagonal. On the other 
hand, for the disordered network the contour lines "circle" around specific links that act as "through-ways" of the network. 
This effect prevents the distance between any two vertices from ever becoming large, that is, of the order of the system size. 
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FIG. 2. Determination of the crossover size n*. (a) 
Semi-log plot of I versus network size for two representative 
values of p and for z = 20. Following Eqs. (0-^), we can de- 
termine n* — apart from a multiplicative constant — from the 
asymptotic slope of £ against Inn. (b) Scaling of n* with p 
for the three values of z discussed in the text. The curves for 
z = 20 and 30 have been shifted up so as to coincide in the 
region where they scale as a power law. Following Eq. (^|) , we 
make a power law fit to n*(p) for p <C 1 and obtain r ~ 2/3. 
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FIG. 3. Data collapse of £(n,p) for z = 10 and different 
values of p and n. (a) Plot of the scaled average distance be- 
tween vertices i/n* versus scaled system size njn* . (b) Same 
data as in (a) but in a semi-log plot. Note the linear behavior 
of the data for n < n* and the logarithmic increase of £ for 
large system sizes. 
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